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We investigate population dynamics in iV-level systems driven beyond the linear regime by a 
strong external field, which couples to the system through an operator with nonzero diagonal 
elements. As concrete example we consider the case of dipolar molecular systems. We identify 
limiting cases of the Hamiltonian leading to wavefunctions that can be written in terms of ordinary 



exponentials, and focus on the limits of slowly and rapidly varying fields of arbitrary strength. For 
rapidly varying fields we prove for arbitrary iV that the population dynamics is independent of the 
sign of the projection of the field onto the dipole coupling. In the opposite limit of slowly varying 
fields the population of the target level is optimized by a dipole resonance condition. As a result 
population transfer is maximized for one sign of the field and suppressed for the other one, so that 
a switch based on flopping the field polarization can be devised. For significant sign dependence 
the resonance linewidth with respect to the field strength is small. In the intermediate regime of 
moderate field variation, the integral of lowest order in the coupling can be rewritten as a sum of 



terms resembling the two limiting cases, plus correction terms for N > 2, so that a less pronounced 
sign-dependence still exists. 
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I. INTRODUCTION 



At] 



When strong few-cycle-, one-cycle- or sub-one cycle pulses [l|-l5|, or arbitrarily shaped 
pulses interact with atomic or molecular many-level systems, significant population transfer 
may occur within a fraction of an optical cycle Concepts based on averaging over 

field oscillations [loj], frequencies and detuning lose their importance, while properties like 
the carrier-envelope phase (CEP) [2| and interactions with permanent or induced dipole 
moments of the system take precedence. 

We have previously investigated population transfer in a dipolar molecular system induced 
by one- and sub-one cycle pulses and found a strong dependence of the dynamics on the 
CEP and in particular on the sign of the projection of the electric field onto the difference of 
permanent dipole moments of the states involved in the reaction, in the following referred to 
as "sign-dependence" of the field S|. A similar sign-dependence for vibrational excitation 
by half-cycle pulses was found by Korolkov et. al. [6j, but not related to permanent dipole 
moments. On the other hand, Doslic et al. ll| and Naundorf et. al. 12| phenomenologically 
discuss a dipole-moment induced tunneling resonance in a 2-level system (2LS) under a dc 
field, but do not go into detail and do not address the case of general pulses. Note that 
earlier Thomas 13] had already given the analytical solution for the population dynamics of 
a 2LS under a constant field. Tatic and Doslic 1^] describe an analogous tunneling process 
in a dipolar molecule driven by a long single field lobe, which they interpret as a distorted 
dc field. In a series of papers, Meath, Power, Brown and coworkers [15l-l20i| consider the 
interaction of laser pulses with dipolar molecules within the rotating wave approximation 



21] , including also the case of a 2LS with a pulse and an added static field |15| . 

In the present paper we take up these points, generalizing the ideas of Tatic and Doslic 
to distort a constant field to lobe-like pulses and putting them on a firm footing. Thomas 



in ref . 



13] did not address the dependence of the dynamics on the field strength and the 



resonance properties of the Rabi-type process. These can be related to the Stark effect 
22[ inducing a change of the energy eigenvalues by coupling to the field, and to WKB- 



like arguments 22|, |23[ suggesting that "resonant" transfer between two levels should be 



maximized at conditions corresponding to degenerate eigenvalues. In this spirit we look 
for conditions for resonances and investigate the influence of the sign of the field; for the 
important point of constructing propagating pulses see below. 
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In the domain of short stron g pu lses traditional tools, like the rotating wave approxima- 
tion [2l| or Floquet theory 24j, |25] , become inapplicable, while methods like semiclassical 



strong field theory 



26 



30t ] are suitable approximations, although only so for large quantum 



numbers. Although exact representations of the population dynamics in strong fields have 



been discussed 



31 



32] , due to their complexity it is hard to gain the physical insight required 



e. g. for constructing a simple field that selectively populates the target level. 

Rather than aiming at exact solutions we search for limiting cases that may lead to 
sufficiently general analytic results for dipole-moment driven dynamics in iV-level systems 
(NLSs) in order to find criteria for effective population transfer. We first note that the 
reason why analytic solutions of the Schrodinger equation do not exist in the general case is 
rooted in the non-commutativity of the two operators in the Hamiltonian corresponding to 
the unperturbed energies and the coupling to the perturbation. This suggests that a way to 
obtain a solution in form of an ordinary exponential, in contrast to the usual time-ordered 
exponential 22], is to search for cases where in a suitable representation of the Hamiltonian 
one of the two matrices can be neglected, and copies of the remaining matrix at different 
times commute with each other. 

We do not further pursue the well-known cases of the weak field limit, which can be 
treated by perturbation theory 22J, and the strong field limit 



26 



3(J. 



Instead we concen- 



trate on two alternative situations, which are characterised by slowly and rapidly varying 
fields of arbitrary strength. In relation to standard frequency-driven conditions both these 
limits would correspond to extreme detuning. The limit of rapidly varying fields can be 
qualitatively analysed for arbitrary iV by transforming to the interaction representation. 
The resulting integrals are elementary and we find that in this limit population transfer in 
NLSs does not depend on the sign of the field. However, for propagating pulses the possible 
population transfer is negligibly small. 

The more interesting case is the one of slowly varying fields, where within certain time- 
intervals the field can be well approximated by a constant value. This case can be solved by 
the adiabatic approximation [33j and diagonalization. We show that in this limit a resonance 
emerges which enables effective population transfer and will be interesting for applications. 
This " dipole-resonance" determines the magnitude and sign of the field in contrast to the 
resonance condition on the frequency in usual spectroscopic transitions under weak fields. 
Furthermore, for 2LSs the reference case of a constant field has an exact analytical solution 
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13| . so that a combination of these results with the present allows a more comprehensive 



understanding of the dipole resonance. 

In the intermediate case of moderate field variation no general analytic treatment is 
possible. In this case we use the contribution of lowest order in the interaction picture 
and show that it can be rewritten as a sum of two terms, each one representing one of 
the above limiting cases of field variation. Thus we obtain the result that a somewhat less 
pronounced dependence on the sign of the field may occur, which becomes manifest only 
after a certain " induction period" . We extend our treatment of this case to 3-level systems 
(3LSs), whenever possible also addressing generalizations to the case N > 3. 

Propagating pulses show the property that in the far-field the time-average of its electric 



field goes to zero 



fly, 
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35]. Such pulses are naturally obtained from vector potentials with 



the property lim t ^. 0O [74(t) — A(— t)] = 0. More pragmatically, "effective half-cycle pulses" 
obeying the zero time-average can be designed as fields with one pronounced lobe balanced 
by long, but weak tails with opposite sign of the field strength or by series of small lobes 
5j. When using single lobes, we anticipate that such small "side" lobes will make negligible 
contributions if their peak field strength is far from the resonance field strength. This 
standard approach of working with single lobe fields is justified a posteriori by simulations 
on model systems [361 ] . 

Our paper is organized in the following way. In section[TT]we present general developments 
of the theory, which lead to the specification of the limiting cases of "slowly varying" and 
"rapidly varying" fields. Section HTT1 is devoted to the case of slow field variation, in which we 
include separate subsections on 2LSs and NLSs. Next section HVl investigates the theory for 
rapidly varying fields, and in the second subsection we combine the results for the limiting 
cases to obtain results for the intermediate case of field variation. In Sec. |V] we give a 
summary of our investigations together with our conclusions. In three appendices we in 
turn discuss spectral properties of propagating pulses, estimate the magnitude of multiple 
integrals of the field, and give details of the derivation of the results for rapidly varying 
fields. 
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II. THEORY: GENERAL DEVELOPMENT 



We consider a sequentially coupled, potentially branched NLS, representing e.g. vibra- 
tional levels of a molecule with a permanent dipole moment. Using the semiclassical dipole 
approximation, the Hamiltonian can be written in algebraic form, and the time- dependent 
Schrodinger equation becomes 



id t c k = [e fe - fJL kk E{t)] ^kiE{t)ci 

l^k 



(1) 



The Cfc, k = 1,...,N, are the time dependent expansion coefficients of the eigenstates 
of the potential with eigenvalues e k . E(t) is the projection of the electric field onto the 
dipole operator with expectation values The diagonal elements /i kk represent the dipole 
moments. 

To proceed we employ a Taylor expansion starting from the diagonal elements of the 
Hamiltonian, and switch to the interaction picture at t = 0. Defining new coefficients 

ft 



we obtain 



<f) k (t) := exp l^i J dt' [e k - /i kk E(t')) J c k (t) 



(2) 



(3) 



where the time-dependent matrix C(t) is given by 

C kl (t) := -fi k iE(t) exp jz jf dt' [Ae kl - A/i kl E(t')} j 
C kk {t) := . 



(4) 



Here we also used 

Ae k i := e k - ei , A/i fc / := ^i kk - fi u . (5) 

Note C k i(t) is non-zero only for coupled pairs k,l. The solution to this problem is the 
time-ordered exponential [22] , 



«(*) = E 



d*iC«(*i) + H) 2 / dhJ^Cki'ih) dt 2 C in (t 2 ) + . . . 
Jo Jo ,, Jo 



■ £(^ ) +4; ) +4? + ---M(o). 



(6) 



In this notation represents a term with an n-fold integral. 



In the following we consider a transition from an initial state i to a final state f, for which 
there exists a unique shortest coupled path with s steps. From eq. (JHJ) we have 



-I 



/ dtiC fl (ti) + H) 2 / dhJ^Cuih) dt 2 C l r 1 (t 2 ) + 
Jo Jo „ Jo 



0i(O) 



= (/« + /f + . . . )0 i (O) , (7) 

where Ig is an r-fold integral corresponding to an r-step process, and hence li = for 
r < s. We define the population Pkif) of a given level k by the square modulus of the 
wavefunction 4>k(t) for this level. Using P;(0) = 1, the population of the target level f then 
becomes 

p f (t):=i^(t)i 2 = X;ii4 n) i 2 +2E Re N n) r4 w) ll • (8) 



Terms of the form (1^ yij? 71 \m ^ n, are interference contributions. 

The properties of the series of time-ordered integrals are well studied [221 ]. yet so far 
no simplifications have been derived for the case of general perturbations E(t), e.g. by 
removing the time ordering or by deriving expressions in terms of elementary functions. 
Note that not even Ig can be treated analytically for general functions E(t). Integration 
of the Schrodinger equation without going to the interaction picture is not any simpler, and 
again leads to time-ordered integrals. Note however, assuming the equation can be solved by 
diagonalization, that the interaction picture does not lead to the same eigenvalues than the 
Schrodinger picture. In order to obtain the correct dynamics, it is essential to diagonalize 
the Hamiltonian matrix of the original Schrodinger equation. 

With this situation in mind we address the question which conditions on the parame- 
ters would lead to an analytic solution. Certainly in the absence of time-ordering eq. (JTj) 
would become an exponential of the integral of the matrix C, and therefore the level pop- 
ulations could be calculated explicitly. Time ordering arises from the noncommutability of 
the Hamiltonian matrices taken at different times. In order to obtain "simple results" we 
have to identify conditions under which these matrices do commute. This leaves us with the 
following cases: 

1. Strong field limit: the energies involved are much smaller than the diagonal contribu- 
tions from the field, e& << ^ikkE{t). The dynamics can be obtained by diagonalization 
of the dipole matrix. 
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2. Weak field limit: The off-diagonal terms in the Hamiltonian matrix are small. This 
case can be treated by perturbation theory [22| . 

3. Slowly varying (adiabatic) field limit: For every t using the expansion of the field 



field E(t) = J2n=o d t E ( t o)( t - t ) n /n\ we assume that \E(t) - E(t )\ « \E(t)\. In 



this case the dynamics can be well approximated by a system of piecewise constant, 
time-independent Hamiltonians. 

4. Rapidly varying field limit: Considering the propagation matrix in the interaction 
picture, we find a formal analytic solution in case of rapidly varying fields. 

In the following we will concentrate on the last two cases. To this end, as a measure of 
the variation of the field with time we introduce Q m { n and f2 max as a characteristic lowest 
and highest Fourier frequency of E(t), suitably determined from the spectrum. Now we can 
distinguish two limiting cases, quantifying in turn the conditions in items 3 and 4 above, 

(a) (slowly varying) min {/M} |Ae H | >> |fi max |, 

(b) (rapidly varying) max W} |Ae w | << |O min |, 

as we show in the following. Here {k, 1} denotes all pairs of levels within the reaction path. 
In terms of frequency, (a) and (b) correspond to the two opposite regimes of large detuning. 

III. SLOWLY VARYING FIELDS 

A. Implications for slowly varying fields 

In order to treat the limit of slow variation, we start by representing the field as an 
expansion around a fixed value to, 



holds, i.e. we deal with "almost constant" fields. In the following we show the implications 
of this relation on the parameters of the field. 




(9) 



n=l 



Due to the low frequency of the field we expect that 




(10) 
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We rewrite the general spectral representation, eq. (1A15I) from appendix |Al with respect 
to inversion symmetry, so that 

i r°° 

E(t) = E(t ) + - / dQS u (Q) sm[Q(t - f„)] + S g (n){cos[Q(t - <„)] - 1} . (11) 

7T Jo 

Here S u (S g ) denotes the spectrum corresponding to functions of odd (even) symmetry, 
and we define the phase, introduced in appendix El as <ft := — ^o- From the Fourier 
representation we explicitly took out the constant term, given by 



7T 



oc 



E(t ) = - I dns g (n) . (12) 



o 



Replacing the sine and cosine functions by their Taylor series leads to 

E(t) = - / dtl {S u (tt){tt(t - 1 ) -tt 3 (t-t ) 3 /3\ + ...} 
n Jo 

+s 9 (n)[i-n 2 (t- to) 2 /2! + ...]} , (13) 

and therefore from eq. and f lT5|) . using E^°\t ) := E(t ) for compact notation, we obtain 
the relation 

\E^(t )\ = - dQS R (ft)Q n , (14) 

Jo 

in which R denotes the type of representation the spectrum belongs to, R = g for n even 
and R = u for n odd. Due to the discussion in Appendix [A] we can translate the notion of 
"almost constant" to the property that the spectrum of E(t) is sharply peaked at a small 
value Oq, even in case we include the switch term. The linewidth of this peak is given by 
2f2o/ct7TA/ln(2), which goes to zero with fio — > 0. Assuming the spectrum is similar to a 



gaussian peaked at Q = 0, from eq. ([H]) we find |£ (n) (t )| ~ (fi /avr v /ln(2)) n |£ , (t )|, which 
gives rise to the upper bound \E^(t )\ < Q,Q\E(t )\. Using this relation in eq. © we obtain 
a geometric series as a majorant for which converges for \Q (t — to)\ < 1, so that we 

obtain 

1^' < ' (15) 

Investigating the dynamics within a given time-interval of length T, e. g. the length of 
the pulse, we deduce from eq. (fT5|) that the condition for slow variation, which is equivalent 
to E(t) w E(t ), becomes 

T « l/Vlo . (16) 
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For the approximation to be true this sets a limit on the integration time, so that to each 
slowly varying field there is a longest time interval beyond which higher derivatives of the 
field become significant. 

Next we derive a necessary condition for eq. f fTUj) . which is often intuitively connected 
with slow variation. We obtain 



\E{t)-E{t ) 



n=l 



(t - t )E(t ) 



n=2 



\E'(t )(t-t )\ , 



(17) 



where we used \Q (t — t )\ « 1. From \E'(to)(t — 1 
that the condition of slow variation is equivalent to 

E(t ) 



T « 



E'(t c 



\E(t)-E(t )\ « \E(t )\ it follows 



(18) 



which is necessary, but not sufficient for eq. flTOjl . 

We pause for a moment to discuss a subtle point about our treatment. In general we 
are interested in results for situations where the field is zero before and after the pulse. 
The present limiting case is unable to deal with a small field strengths (at the tails of 
a pulse), because the derivatives will eventually exceed the field strength even for slow 
variation. However, the case can be extended to complete pulses if the tails do not contribute 
significantly to the population dynamics. In numerical studies [36[ we indeed find this to 
be the case for the 2LS, and for N > 2 with some restrictions due to alternative transfer 
pathways. Furthermore, by studying the adiabatic approximation of a general 2LS and non- 
adiabatic coupling (NAC) to first order we find that the slope of E(t) in the tails of the 
pulse, which gives rise to the NAC, has hardly any effect during the rise of the field strength, 
but is responsible that during its decrease the population remains in the target level [37]. 



B. The two-level system 

We now address the simple case of a 2LS, where the concept of a "dipole resonance" is 
particularly clear, and where the reference analytical solution is available for a 2LS under 
a constant field as the ultimate limit of slow variation. This section is hence connected to 
the question if, and how, a resonant constant field may be deformed into a pulse so that 
maximum population transfer is maintained. 
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1. Dipole resonance in the two-level system 



In the 2LS we have only one possible step from level 1 to 2 (which in our notation now 
become i and f). In order to use the results of this section for N > 2 discussed below, we 
consider the integral M% corresponding to the step i to f as a possible innermost integral of 
an arbitrary term in eq. ([7]) that is part of a transfer pathway in an iV-level system. 

By using the approximation E(t) ~ E(t Q )[l + (t — t )Q + 0(Qq)] we obtain 

/ dtE(t) = E(T) dt[l + (t-T)Qo] =TE(T)[1 + — -] mTE(T) (19) 
Jo Jo 2 

for the integral in the exponent of Mi . This is similar to the adiabatic approximation 22|, |33) , 
where the field is treated as constant when performing the diagonalization and integration 
of the Hamiltonian, with subsequent substitution of E(t) for E in the result. We get 

Mxih) « m.iEitx) [ dt exp {i [Ae {ll - A/i hi E(t )} t } 
Jo 

„ HuEjh) [exp {i [Ae hi - Afi^Ejh)] h} - 1] 

This equation clearly shows a resonance whenever E(t) = Ae^i/ A/iy =: Aq, which we term 
dipole resonance. We note that the resonance selects one sign of the field. 

2. The two-level system in a constant field 

The resonant behaviour seen in eq. ( 12"U|) is equally obtained for the 2LS under a constant 
field, as expected. As noted by Thomas [13], the analytical solution for the evolution of 
population in the 2LS under a constant polarized field with (effective) field strength E is 
governed by Rabi-dynamics. It is convenient to use reduced parameters, slightly different 
than the ones used in which are defined as 

d := A£fi ~ ^ E , 6:=^t. (21) 
2fiaE '2 v ; 

We call d the dipolar detuning, zero at dipole resonance, and 6 is a generalized time. Using 

the initial condition P\(t = 0) = 1, we have 

PliB) = sin 2 (eVTTa?) . (22) 

We observe Rabi-like behaviour with 9, where the dipole resonance is reflected by the pref- 
actor containing the detuning in form of a Lorenzian. However, as shown in Fig. 1, variation 
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of the "physical" tuning parameter E does not result in a Lorenzian. We plot the maximal 
final population P f max (d), found by varying only 9, against E and d in Fig. [TJ Note that 
d = represents the global maximum of P f ma:r and in the limit of infinite field strength we 
only get P f maa; = ^/ (^u 2 . _|_ A/x|/4), independent of the sign of the field. Using the resonance 
field with inverted sign, E = — Ae^/Afi^, we obtain a maximal population in level f of 
pmax _ ^2y^2 _|_ A/ig). Immediately we conclude that significant sign dependence occurs 
only for yU fi << A/i fi . 

Next we use this criterion to derive a condition of sign- dependence for slowly varying 
fields. We map the slowly varying field onto a constant field by averaging E(t) and then 
combine this relation with the one for slow variation in eq. f)16p . Assuming we constructed 
E(t) such that the average field strength is on resonance we obtain from eq. f[2"2"j) that the first 
instance of population inversion is given by the Rabi time t = tt/ [irAq. Inserting into eq. ffT6j) 
we conclude Q << (/ifi/7rA/ig)Aefi. Hence we have two possibilities to fulfil this condition: 
either we make Q very small or we search for systems with /ig >> A/ig. In the latter case 
Q need not be too small but there is no sign- dependence. We note that only in the former 
case do we get sign-dependence for which we find Q << Aeg as a necessary condition, which 
relates system parameters to field parameters. Note that, supported by fig. [TJ a significant 
sign dependence implies a significantly narrow resonance peak as a function of E. 



C. iV-level systems 

In order to illustrate a general system with s > 2, we use the example of a 2-step process 
i— >b— >i from an initial state i via an intermediate state b to the final state f in a 3-level 
system. 

(2) 

We start by analyzing 7 fi , which is the simplest term that already reflects the added 
complexity of more than two levels. Jg is the term of smallest order in the coupling that 
contributes to the population of the target level. Explicitly we have 



4 2) = H) Wm / dti / 1 dhEit^Eih) 
Jo Jo 

xexpji Aefbti + Ae bi t 2 - J dt' ' A^E(t') - J dt' Afi hi E(t') j 



Consistent with the condition of slow variation of the field we approximate this integral in 
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the following way, similar to the procedure in Sec. IHIB 1\ 

4 2) ~ H) Ww f dti / 1 dt 2 E(ti)E(t 2 ) 
Jo Jo 

x exp {«[A efb - A// fb S(t 1 )]t 1 + [Ae bi - A/i bi E(t 2 )]t 2 } . (24) 

From our discussion above we know that the tails of the pulse can be ignored and therefore 
we use E(Q) = as the lower boundary of each integration. This gives 

r(2) ~ ,, ,, F(f )2 exp(t[Ae fi -A/^E(t)]t) 

[Ae fi - A/Xfl£?(t)J[Ae b i - A/j b i-E(i)J 

where Aeg, A/ig are differences between values at i and f, and Ae b i, A/i b ; describe a transition 
from i to b. We assume for the moment that E(t) ^ Ae b ;/A/i bi holds, i.e. no resonance 
occurs for intermediate steps. Eq. (1251) then clearly shows a resonance at E(t) = Aq = 
Aefi/A/ifi, in which case the population of the final state becomes 

P&\E=Ao(t) = filfilA^/lAeu - Afi hi A ] 2 . (26) 

We call this an overall resonance because it establishes a resonance between the initial and 
the final state across all intermediate steps. 

We return to the case of simultaneous resonance of any number of intermediate steps 
(but not all of them). In our example this reduces to the single step i— >-b, so that we obtain 
the condition E(t) = Ae bi /AyU bi . Eq. ( 1M|) then becomes 

4 2) w H)VfbyUbi£(t) 2 J dhti exp |i Aefcti - J dt'A^E{t') J . (27) 

The remaining integral contributes significantly only if an additional resonance condition is 
fulfilled - in our example the resonance of step b— >f. Apart from the case that the field was 
shaped to attain the correct resonance condition of each step just at the time transfer occurs 
there, this is only possible for systems where all steps have the same A (" dipole-harmonic" 
systems). For this case no distinction with respect to N is necessary. In dipole-harmonic 
systems transfer to any "final" state will not be highly selective as all intermediate levels 
are on resonance. 

Our analysis for a two-step process has a natural generalization to s > 2. Dynamics via 
an overall resonance (direct coupling of i to f even if /xg = 0) is direct and hence different 
from step-wise population transfer. From the results obtained above for the 2LS and the 
3LS we conclude that population transfer is maximized when an "overall" dipole resonance 
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condition is met. The resulting resonance field is determined by a ratio of the average of 
energy differences and the average of dipole moments along the shortest path P between the 
initial and the final state, 

A„=fs^ = ^. (28) 

We note the interesting detail that within this approximation the parameters of the interme- 
diate levels have no influence on A Q (which reflects overall resonance). However, Pf depends 
on the properties of the intermediate level by eq. (I26j) . 

IV. RAPIDLY VARYING FIELDS AND THE INTERMEDIATE CASE 
A. Rapidly varying fields 

The limit of rapid field variation is defined by the condition max^/} |Ae^| << |f2 max |. 
We show in Appendix O that this limit is applicable whenever the majority of the Fourier 
spectrum of E(t) is located well beyond the largest energy difference from the initial level to 
any level along the transfer path. An example would be a single field lobe with sufficiently 
short duration (there is no immediate restriction on the number of lobes). The n levels along 
the transfer path are denoted {/ , h, ■ ■ ■ , l n } with l Q = i and l n = f . We investigate the first 
step from i to l\ and begin with defining the auxiliary functions 

a(x) := exp(zAe/ li x) , b(x) := E(x) exp ^— zA/i/^ J dx'E(x')\ 

B(x) := f dx'b(x') (29) 
Jo 

Due to the oscillating kernel, repeated integrals of b could be zero at isolated points. In 
order to keep our results general, we consider the repeated integral of order k = n nonzero, 
while all integrals of lower order with k < are zero. Using uq partial integrations, in 
Appendix [C] we show the validity of the approximation 

Mi(ti) = i/i^i / dxa(x)b(x) 



o 



(-l) n °-V w ^ no ^{x)I B {n ,x)\^ - J dxa^\x)I B (n ,x) 
z/i ili (-Ae Zli )" - 1 a(t 1 )/ B (n 0) ti) . (30) 
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in the present limit. In eq. ( 130]) we use the definitions Is(k, t) := L dxlsik — l, x); Ib(1, t) = 
Bit) as the fc-th integral of Bit), and denotes the k-th derivative of a. The properties 
of fc-fold iterated integrals of the field, Ik(E,t), are discussed in appendix iBl 

Only two cases are relevant when considering pulses: the integral over the pulse may be 
zero (propagating pulse) or nonzero, leading to n = 2 or n = 1, respectively. 



1. n = 1 



Explicitly inserting a and b into Mi and using eq. ( 1C7I) . we obtain 



Mi(ti 



— exp(zAe w ti) <H - exp 



iA/i^i / dt'P(t' 



(31) 



(2) 

We use this expression in the kernel of the integral l l ■ £ that represents the 2-step process 
(from i to Z 2 ) and apply eq. f )30|) once more to obtain 



M 2 (t 2 



^W e xp[i(Ae Wl +Ae w )t 2 ] 



1 



l-exp[-iAn hh [ dt'E(t')} 
Jo 

rti 

1 - ex P H(A/i Wl + Afi hi ) / dt'E(t')} 

Jo 



(32) 



Proceeding in this way we note that contains a sum of terms with the same functional 
form as Mi, which however depend on Ae^i and A/x/.i and correspond to a subpath of the 
reaction path from i to lj, with j < k < n; for details see Appendix O 

The calculation of the target state population Pf involves taking the square modulus of 
a sum of terms consisting of a real factor multiplied by a product of the complex quantities 
Pk{x) and f/3a, which are defined in eq JCll It follows that Pf is written as a sum of factors 

{Pkfa^)*Plfnu + PkfaftiPlf K u)* = 2 < COS ^ Ae h+ih x 



COS 



+ COS 



fc-1 



J 

'fc-1 



r dfE(t') 


— cos 


° 





y^Ae lj+ll .x + w w [ dt'E(t') 
3=1 Jo 



IV i „ x 

V, Aei^x - (u^ - ty Ki ,) / dt'E(t') 



(33) 



multiplied by real numbers. Without loss of generality we used n > k > I. 
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Regarding eq. (|33|) . a few comments are in order. In case of a single initial state in a 
system with a unique path to the final state the sum over energy differences does not appear 
in eq. ( I33p . because direct contributions to the wavefunction and interference terms (path i 
to f augmented by loop paths from f back to f) all contain the same phase factor of energy 
differences. We arrive at the important conclusion that in this case the population of the 
final state does not depend on the sign of the external perturbation E. Note this result 
follows without any assumption about the number of field lobes. 

If the initial population resides in more than one state, more than one path to the target 
level contributes to the dynamics, each one with its own energy difference and corresponding 
phase factor. In general, due to interference effects the dynamics will then depend on the 
sign of the field. This sign dependence however vanishes asymptotically in the strong field 
limit. 



2. n > 1 

This case occurs whenever the field is supplied as a propagating pulse. It implies that the 
time-integral over the field is zero 



35|, resulting in 5(0) = 0. This is a common situation in 
experiments which we therefore discuss here separately. In the general case we arrive at the 
estimate 



no-1 



= 0[ l ™ ) . (34) 



Equipped with this relation of general order we investigate population transfer at the end of 
a propagating pulse, Pf(t — > oo). We can apply the approximation in eq. ( 1301) with no = 1 
for all integrations but the one corresponding to the last step, which contains t — > oo as an 
upper limit. The reason lies in the fact that B(x) becomes zero at x — > oo (for all times after 
the pulse has passed). Hence we have to consider n = 2 only in the integral representing 
the last step to level f. The same holds true for all possible interference terms. Thus from 
eq. (IC4I) we expect the maximum final population to be of the order (Aen^/fio) 2 << 1 
smaller than in case no = 1. Due to the fact that this holds for arbitrary s, we conclude that 
a short, rapidly varying propagating pulse hardly transfers any population. It is important 
to recall that rapidly varying in our definition refers to the total field E(t), independent 
of the envelope, and hence to a field of very high frequency. Like for n Q = 1 we find that 
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population transfer (regardless of its small magnitude) is independent of the sign of the field 
if we start from a single initial state. 



B. Intermediate case 



1. 2-level systems 



If |Aefl| is comparable to both f2 m i n and n max , we are in an "intermediate" regime con- 
cerning the two limits considered above. No analytic treatment is available, and we have to 
use other tools to analyse the integrals of the time-ordered series. We first note that due 
to the same functions E(t) occuring in the exponent and in the factor multiplied with the 
exponential in the kernel of the integral, it is useful to use the identity 



E(t) = -(g' kl (t)-Ae kl )/A l 2 



ki 



with the time- derivative of the phase defined by 



g' kl (t) = Ae kl - Afi kl E(t) 



(35) 



(36) 



The indices k, I are arbitrary and correspond to a given step. Using this relation in the 
integral representing the first step together with substitution leads to 



Afi 



fi 



cxp 



Ae fi t 



/ Aeg I dt exp < / 



AfiQ / dt'E(t' 
Jo 

/•to 

Ae fi t -A/i fi / dt'E{t' 
Jo 



1 



(37) 



The first term is equivalent to the expressions obtained in case (a). The second term can 
neither be calculated analytically nor be approximated as slowly varying in the whole domain 
of integration. However, from the discussion in section IHIB II it is clear that for sufficiently 
short integration time we could apply the slowly varying approximation to the integrand. 
This leads in a natural way to the idea of partitioning the domain of integration into intervals. 
It only remains to find the intervals that show resonance. In the optimal case only one 
significant contribution remains and we can indeed replace the second term in eq. ( 1371) by 
its slowly varying approximation. 

The integrand of the second term in eq. (1371) cannot be treated by the saddle point 



approximation (SPA) 



because H — > need not be true and in addition we integrate over 
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A > A > 


A 


> > A > 


-A 




condition: 


or 




or 




\A\ > \A \ 




A < A < 


A 


< < A < 


-A 




solution: 




X 37T 

d 2 - 27J 


h,4 = ^arcsin(^) 



TABLE I: Possible solutions of eq. (|38p in [0, 2ir] corresponding to minima of g for E(t) = Asinut 

a finite time-domain which results in contributions from the contour near to = and to — t\. 
Elsewhere we will discuss a method of partitioning the domain of integration 



To find resonance points 5, we search for a minimal first derivative of the modulus of the 
phase in the exponent. Around these points the variation of the phase factor in the integrand 
is slow and the resulting integral becomes large. Note this is a more general criterion than 
in the SPA, including all possible saddle points. 

In the following we abbreviate the phase difference between initial and final state with 
g(t). In order to use differentiation we note the bijective mapping of \g'\ to g' 2 . This leads 
to 

d t [g'(t)] 2 \ t=s = 2g'(S)d t g'(t) \ t=8 = , (38) 
and demanding a positive second derivative, 



d 2 t W(t)f 



\t=s 



2 [d t g'(t)f 



\t=s 



2g'(5)dfg'(t)\ t=s >0 



(39) 



we obtain the centers 5 of the resonance intervals with maximum population transfer. If 
\g'{t)\ > for all times, then 5 lies at an extremum of the field. Otherwise g(t) becomes sta- 
tionary, g'{8) = 0, so that E(S) = Aeg/A/ifi =: A , which we denote true "dipole resonance". 
The latter case corresponds to the usual condition of the SPA. The above equations show 
that for oscillatory fields population transfer can only be large for one given sign within the 
period. The opposite sign corresponds to a maximum of g' 2 and yields hardly any transfer. 
In passing we note that the optimal case of a dipole resonance can only occur if g(t) is not 
strictly monotonic. 

In general, we may obtain more than one solution 5. As an example, for E(t) = A sin cut 
the solutions of eq. ( 1381) that lead to maximum population transfer are given in Table [H If 
\g'{t) \ remains sufficiently small between two solutions, the corresponding resonance intervals 
merge, and the dipole resonance condition is maintained for a particularly long time. If E(t) 



17 



consists of a single half-cycle lobe, its optimal amplitude A should therefore be somewhat 
larger than Aq. 

From eq. (j37p it is apparent that the resonance condition has to be maintained for a time 
longer than l/Aeg, or equivalently 1/t < Aeg, in order that the second term dominates the 
sign-independent first term. Combining this result with the condition of a slowly varying 
field, Qt << 1, (at least valid around the interval of resonance) gives a condition relating 
properties of the field to the system parameters, namely mini^} |Aen| >> |^max|- Note this 
relation is identical with the one presented in section IIII B 21 Furthermore, it is opposite to 
the one for fast variation in section IIVI The sign dependence gradually disappears as the 
pulses become shorter. For very short pulses this case goes over to the limiting case (a). 



2. N -level systems 

For N > 2, new aspects arise due to the fact that the single steps along the transfer path 
can have different resonance amplitudes. This necessitates the definition of what we call 
diagonal detuning. For a given (sub)path the diagonal detuning d p can be defined relative 
to a reference transfer path from initial to final state by comparing the phases, 

d p {t) := g' p {t) - g'S) ■ (40) 

We realize that the diagonal detuning is time-dependent for general fields. At a true overall 
resonance we get d p (S) = g' p (S) as the phase difference along the path p. Note in case of a 
dipole-harmonic system d p (5) = for every subpath of the transfer path. 

We again demonstrate the dynamics by discussing a two-step process. Similar to the 
treatment of the case iV = 2 we use eq. (|35|) and substitution to split each of the two 
integrals in M 2 into a sum of 2 terms. We first substitute for E in the inner integral , 

GO = H)2 /^bi /Vi5(f0expfoftb(fi)] 

^Mbi Jo 

x |(exp [ighi(ti)) - 1) - iAe hi [ dt 2 exp [igu(t 2 )} 

1 •\2A i fb/ i bi 







A/ibi 



/ &t\E{t]) exp [igaih)] (1 - exp [-i# b i(*i)D - 
Jo 



- iAe hi dti£(ti) exp [ig^h)} dt 2 exp [igu(h)] j . (41) 
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It is interesting to note that population transfer shows not only the expected resonance from 
the initial to the final level, indicated by g&, but also a concurring one to the intermediate 
level. The latter resonance has the effect that in case of a harmonic system mainly the third 
term contributes and no resonance from the initial to the final level occurs. The third term 
is a correction term with no further possibility of simplification, which describes stepwise 
excitation. Using eq. (|35|) with indices f and i for replacing E(ti) in the first two terms, and 
after writing the phases out explicitly we finally obtain 



A^Afl 



- jexp (i[Ae & t - J dt'E{t')\\ ~ 1 
- zAe fi J dti exp (i[Ae^t x - A// fl J dt'E(t')]\ 
+ zA/ifi jf dti^(ti) exp ^[Aef^ti - A/i tll dt'E(t')} 

x 1+iA^iJ dt 2 exp(i[Ae hi t 2 -Ai Mli J^ dt'E(t')]^ J. (42) 

Apparently, an interpretation similar to the case N = 2 above can be given. The first 
term does not induce any dependence on the sign of the field. The second term gives rise to 
the overall resonance, with conclusions equivalent to iV = 2, if the appropriate differences 
of energies and moments are used. The last term is a correction term and consists of 
two iterated integrals, each one with a resonance field of the respective step. In case of a 
harmonic system this term will contribute along with the second term. In case of systems 
with small but nonzero diagonal detuning a complicated temporal behaviour with beatings 
is possible. In case of large diagonal detuning along the path the third term will be negligible 
and the dynamics should be similar to an effective 2LS. Generalizing to more than 2 steps, 
the number of correction terms will increase. In case of alternate paths only such correction 
terms can contribute to quantum interference effects. We conclude that in systems with 
large diagonal detuning fairly symmetric switching of population by an appropriate choice 
of the field strength is possible because no interference can occur. 

Together with the results obtained for the 2-level system in the intermediate regime we 
conclude that in systems with sufficient diagonal detunings the relation 

remains valid for fields with moderate variation. The quality of the approximation of inde- 
pendence of Aq on the details of the path increases with increasing diagonal detuning. 
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In all situations discussed up to now we have assumed that each A/z relevant for A is 
nonzero. Considering for the moment the case A/Zfi = 0, we see that Aq will tend to infinity. 
This expresses the fact that in this case a resonance cannot occur, and population transfer 
can only proceed through other mechanisms. 

Elsewhere we will present numerical simulations of population transfer in 2LS and 3LS for 



pulses and constant field [36|, which clearly show the qualitative features discussed above. 



V. SUMMARY AND CONCLUSIONS 



Using the time-dependent Schrodinger equation, we analyse the population dynamics in 
iV-level systems induced by strong and short field pulses. The operator of the system, which 
couples to the external field, is assumed to contain different nonzero diagonal elements. In 
the present article we investigate as an example a system with permanent dipole moments 
coupled semiclassiclally to an electric field via the dipole operator. 

We search for possible limiting cases of the Hamiltonian, for which the solution of the 
Schrodinger equation can be well approximated by ordinary iterated integrals without time- 
ordering. Using respectively the original Schrodinger picture and the interaction picture, we 
identify two new regimes relevant for strong fields, namely the two opposite limits of large 
frequency detuning. 

In the limit of rapid variation of the field, compared to a characteristic energy difference of 
the system, we prove for arbitrary N that the population dynamics is independent of the sign 
of E(t). Furthermore, for propagating pulses no significant population can be transferred 
by the pulse. 

In the opposite limit of slow field variation, population dynamics is determined by a 
resonance originating from the diagonal matrix elements of the observable contained in the 
coupling. This resonance selects the amplitude and phase of the field. This is in sharp con- 
trast to the well-known resonance in the perturbative regime that determines the frequency 
of the field as a difference of eigenvalues for a transition. As a result population transfer 
crucially depends on the sign of the projection of E(t) onto the appropriate difference of per- 
manent dipole moments whenever the resonance linewidth with respect to the field strength 
is sufficiently small. 

We also give a qualitative discussion of the intermediate case of moderate variation of 
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the field, by analysing the contribution of lowest order to the population in the final level. 
For N = 2 we rewrite the integral to a form consisting of a sign-independent and a sign- 
dependent term, each corresponding to one of the limiting cases discussed above. A signif- 
icant sign dependence occurs only for field pulses with a duration longer than the inverse 
level spacing of the initial and final level. In case of N > 2, added correction terms must 
be considered. The correction terms usually show no resonance but their number increases 
with increasing number of levels, weakening the effect of the resonance. 

The pronounced sign-dependence of the population dynamics for slowly varying fields 
would allow for ready control of branching between target states with different polarity and 
thus suggests an application of dipole-resonant population transfer in a setup acting as a 
molecular switch. For few-cycle pulses, the behaviour of the intermediate case will also 
become manifest as a dependence of population transfer on the carrier envelope phase. 

The nature of the coupling and the origin of the diagonal terms is not relevant to our 
analysis, and hence our findings should also hold in various problems of optical and magnetic 
population dynamics. Even the condition of existence of diagonal elements of the operator 
coupling to the field can be relaxed by considering induced quantities within an effective 
Hamiltonian of dressed states. This might explain phenomena like the asymetric escape of 
electrons as a non-linear field effect in atoms. 
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Appendix A: Spectral properties of propagating pulses 

An experimentally admissible propagating pulse is formed by shaping the vector potential 
A(t) under the restriction that the difference of A at the beginning and the end of the pulse 
is zero. We assume a gaussian vector potential and write A(t) for the projection of the 
vector potential onto the current density to obtain 



Am(t) 




21 



A(t) 



sin(f2 £ + 0) 



m(t) = exp 



(Al) 
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Here A denotes an amplitude and Qq the frequency. The half-width at half maximum of 
m(t) is given by a/ ^/ln(2). This keeps the discussion fairly general because an arbitrary 
smooth envelope could be well approximated by a finite sum of gaussians with appropriate 
parameters. Using E(t) := —-dtA(t) as the corresponding projection of the field onto the 
dipole moment we obtain 

2(7 — / 



E(t) = Am(t) 



cos(f2o^ + 



QnCT 2 



sin(f2 t + 



(A2) 



We adjusted a to accommodate a optical cycles in the full width at half maximum of the 
gaussian envelope, i. e., a = an^in 2/fi (see also ref. Q). 



Starting from the vector potential in form of a gaussian pulse, in eq. flA2j) we obtain two 
terms for the electrical field. The second term is usually denoted the switch term. We will 
call the first term the principal term. We first derive the spectral representation for the 
principal term, which can be rewritten as 

t 2 



E P (i) = A exp 



a" 



cos(f2 t + 



(A3) 



by changing the origin of the time axis to t max and redefining at the same instance. 
Furthermore, we use flo > without loss of generality. The spectrum can be calculated 
analytically from its definition, 



S P (tt) :-- 



dtE P (t) exp(-int) 



(A4) 



Act— { cos(0) 



exp 



<j- 



4 



— (n - n y + exp -— (n + a 



<j- 



+i sin(0) 



exp 



~(n-no) 



exp 



Aa—^- [Sp t i(Q) cos(0) + iSppity sin 



(A5) 
(A6) 



Clearly the spectrum consists of two parts, which are even (Sp^(Q)) and odd (Sp^ity) with 
respect to inversion symmetry. 

In order to find the spectral representation for the switch term, we again change the 
origin of the time axis and redefine <fi to obtain 



E 3 (t) 



2At 



exp 



t l 



a- 



sin(f2ot + 



(A7) 



Note that the switch term contains a sine, but due to its prefactor belongs to the same 
representation with respect to inversion symmetry as the principal term. The spectrum can 
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be calculated in a similar way as for the principal term, 

poo 

S s (tt) := dtE s {t)exp{-ittt) 



—A- 



2Qn 



(o - o ) exp ( -^-(o - o ) 2 ) - (o + o ) exp ( -^(o + o ) 2 



T 



+i sin(0) 



cos(0) 



(O - O ) exp ( - — (O - O f J + (O + O ) exp ( - — (O + O 



== ~^£r cos (^) + ^(O) sin(</>)] . 



(A8) 



With respect to inversion symmetry, the spectrum again contains an even part (Ss,i(ty) 
and an odd one (65,2(0)). 

Using eq. ( 1A6|) and (IA8[) we find the spectral representation of the field by inverse Fourier 
transformation, 



E(t) := — / dO[Sp(0) + Ss(0)]exp(iOi) 

5s,i(0) 



(A9) 



.4 



dO 



- crSp^O) 



dO 



65,2(0) 
O 

crSp^O) 



O 

sin(0) sin(Ot) 



cos(0) cos(Ot) 



On 



cos(Ot + 



The approximation in the last line uses S_ >av := (Si + S z)/2 and rests on the fact that 
the even and odd spectral contributions from the principal and switch terms form peaks 
that have similar shape near On. Invoking the relation between frequency and energy, this 
approximation can be rewritten as 



E(t) 



Aa 

27^ Jo 



dO 



Sp,cm;(0) — 



cos(Ot + 



(A10) 



an\/\n2_ 

The spectrum of the principal term shows a symmetric peak at On, whereas the spectrum 
of the switch term changes sign at O . It possesses a negative and a positive peak with 
positions 

Os jPe afc = O ± , (All) 



obtained as the local extrema of Ss,au(0). 

From the shape of the spectra and eq. (lAlOj) it becomes clear that the resulting pulse has 
its peak slightly below O and extends farther towards smaller values of O than larger ones, 
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due to cancellation of the principal and switch term spectra. The cancelation introduces a 
cutoff at 

fi c ^ fi + — = fi (i + —^=] . (A12) 

We estimate the position of the peak of the total spectrum and obtain the maximum at 



"'^l^rfT^I (A13) 

(A14) 

The approximation is valid whenever 6a 4 /fig >> 1 holds and fip is very close to fio in this 
case. 

Apparently in eq. (lAlOj) only a cosine turns up because in case = the vector potential 
is an odd function of t and the field will be of even symmetry. The most general case contains 
terms of both characters and can be represented by 

1 f°° 

E(t) = - dfi [S U (Q) sin(fit + 0) + S g (Q) cos(fit + 0)] . (A15) 
Jo 

In appendix [B] we use this form to discuss multiple integrals of the field pulse. We note that 
we could remove the phase by taking = fito, obtaining the usual Fourier composition. 



Appendix B: Multiple integrals of oscillatory functions 

For the proof in section [IV] we need to estimate the size of repeated integrals of the 
function b(x) := i^i^Eix) exp (— iA/i^; J* dx'E(x')) . Note in case of significant variation of 
E(t) the phase factor is close to 1 as its exponent is small. Due to \b(x)\ < \fii 1 i\\E(x)\ an 
investigation of respective integrals of E(x) can be used to find an upper bound. 

We write the spectrum of E(x) as a sum of an odd and even part, S(VL) := S U (Q) + S g (Q), 
and consider the single integration problem, 

pt -| /*oo t*t 

I dxE(x) = - dO / dx[S u (Q)sm(nx + (j)) + S g (n)cos(nx + (j))} 
o ^ Jo Jo 

Here denotes a possible phase and we used the well-known fact that the integral of sin (fix + 
0) cannot be larger than the contribution of half a period. However, already at this point it 
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is clear that this estimate is applicable only if lim^o ^X^) goes as Q n with n > 1. We will 
generalize and refine this estimate below. 
We are interested in the repeated integral 



Pt PXn-1 

I n (E,t) : = / ... / dxi . . . dx n E(x n ) 
Jo Jo 



■y POO pt PXn-l 

- I dQ ... dx 1 ...dx n [S u (Q)sm(Qx n + (j)) + S g (Q)cos(Qx n + (f))] (B2) 

^ Jo Jo Jo 



Using the following exact expressions, 



/ drr sin(f2:r + 0) = — [cos(f2t + 0) — cos(0)]/Q 
Jo 

ft 



/ 

Jo 



dx cos(Qx + 0) = [sm(Qt + 0) - sin(0)]/fi , (B3) 

'o 

we can perform the integrations over time. Iterating we obtain the "boundary term" from 
the lower limits of the n-fold integral, given by 

B n = B^r/n + SiO) -B(n,0)/fi n (B4) 
S(n) := \s u (Q) -S g (Q)} (B5) 



. lWnl fsin(0) — cos(0)l T n even 
B(n,0) := (-l)*^- 1 )/ 2 ! { [ Vr; yvn 

[cos(0) sin(0)] T n odd 



(B6) 



when defining r := Qt and the spectral and boundary vectors S(f2) and B(n, 0) respectively. 
Using the boundary contributions we obtain the following general expression 

I n (E,t) = - [ dQS u (n)F u (Q,t) - S g (Q)F g (n,t) (B7) 



F u (n,n,t) :-- 



<_^irU[(n+l)/2] I Sill «( fit + <P) n eVen 

I cos n (f2t + 0) n odd 



,, . cos n (f2t + 0) n even 
F g (n,Q,t) := (_l)-*[(«+2)/2] <J "V W _ (Bg) 

sin n (f2t + 0) n odd 

where the odd and even functions F u and F g have been defined in terms of what we de- 
note modified sine- and cosine-functions. The first few of these functions and the general 
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expressions are 



sini(r + 0) = [sin(r + 0) — sin(0)] /£l 

cosi(r + 0) = [cos(r + 0) — cos(0)] /Q 

sin 2 (r + 0) = [sin(r + 0) - [sin(0) + cos(0)r]] /Q 2 

cos2(r + 0) = [cos(r + 0) — [cos(0) — sin(0)r]] /Vt 2 



sin n (r + 



cos n r + 



i _ - 93 



j=0 

n.-l 



j=0 



T=0 



T=0. 



sin(r + 0)/n n 



cos(r + <p)/tt n . 



(B9) 



The modified function of order n is constructed by removing the first n terms of the power 
series expansion of the respective sine or cosine function. 

In case r > 1 we find an upper bound for these function as follows, 



sin„(r + 0)| 



sin(r + 



sin(0) + cos(0)r + sin(0)- 



< t r 



i_ ^ i 

An + 1^ A\( r \ 



(r) 



j=0 



j!(r)»-i 



(BIO) 



Within this approximation, the same upper bound applies to cos n (r+0), and for the absolute 
value of both modified functions of order n we obtain the upper bound 



f n 

;i+e)-, 

T 



(Bll) 



by replacing each inverse power of r by 1/r and using the monotonic behaviour of 

with n. The limit r — > oo of the modified functions is smaller but shows the same power in 

t, 

lim | sin n (r + 0)| < t n O([(n - l)!^ 1 ) . (B12) 

T— >00 

For r — > the right side of eq. fIBlOj) diverges like the n-th power whereas the modified 



26 



functions stay finite, namely 



int[(n+l)/2] 



lim sin n (r + 



lim cos n (r + 



cos(0)r + 0{t 2 ) n even 



n+l 

( _ 1)l nt[(n + l)/2] . 

i — - sm(0Jr + U(t ) n odd 



n+l 



-l)int(n/2) t n J s j n (0) + Lll 

nl \ COS(0) + 

.iynt[(n+l)/2] t n J cos (0) + (-1)'"^^ + ^2) „ 

sin(0) 



even 



n! 



C_ 1 - ) int[(n+2)/2] „ 

^ cos(0)r + C(r 2 ) n odd 



n+l 



(B13) 



We note that terms of equal order of the series for the modified sine and cosine function 
differ only by a sign and replacing sin(0) with cos(</>). Therefore we define auxilary functions 



iS n (r),C n (r) to write 



sin n (r + 



^ -^ynt(n/2)-j-n 



+ 



nl 

_-^ynt[(n+l)/2]j.n 

(n + l)! 



sin(0) ( 1 — 



T 



(n + l)(n + 2) 



+ . . 



cos(0) I r 



/ i \int(n/2)j.n 

-sin(0)C„(r) 



(n + l)(n + 2)(n + 3) 

_^-jini[(n+l)/2]^n 



cos(0)5 n (r) . (B14) 



nl (n + l)\ 

in case of even n, with similar expressions for odd n and for cos n (r + 0). From the properties 
of alternating series it is clear that for r < 1 we have |«S n (r)| < 1, |C„(r)| < 1. This leads 
again to a common upper bound 

(n + 2)t r - 



(n + l) 



for t < 1 . 



(B15) 



for both modified functions. 

We now collect the results from eq. (IB 1 1 [) and eq. f)B15j) to estimate the size of I n (E,t). 
We split the integral over Q at r = 1 which yields 



1 4. CM I < 



f 



n-l 



71 



{n + 2) 



dr^r/t) + (1 + e) / dr 
(n + l)! Jo ii r 



S(r/t) 



< 



71 



n + 2 



+ 1 + e 



(B16) 
(B17) 



(n + l)! 

From the definition of the spectrum [see eq. (lAlOp ] and E(t) of a pulse in eq. (1A2f) we get 
-E(O) = Acos(0) < ^4 and the estimate becomes 

n + 2 



\I n (E,t)\<At n 



(n + l)! 



+ 1 + e 



(B18) 
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Appendix C: Rapidly varying fields 



First we assume the approximation in eq. (130|) valid and work out the general structure 
of the terms in 0f by abbreviating the following quantities occurring in M n , 

/ n-l \ n-1 

p n (x) := exp liJ2 Ae i j+ iij x II^+i^ 



j=0 / j=Q 

w Pa := A Vi j+ iij 



-i 



(x) := 1 — exp 



-iwpa [ dt'E{t') 
Jo 



(CI) 



Now we can write the first few M n as 
M 1 = pj'~ 
M 2 = p 2 
M 3 = p 3 



fio 




Wxo 




( /21 


fx \ 


\W 2 lW W 


W20W10J 


f'32 


(--- 




\W 2 l w 20 



/31 /30 



(C2) 



W 31 W 21 Wio W 30 W 20 Wi _ 

It is apparent that there are exactly n terms in each M n . The increase of the number of 
terms with increasing subscript comes from the two parts contained in /, namely 1 and the 
exponential, so that each term in M n can be interpreted as a transition from some level 
lj,j < n to l n . We do not pursue the details of the calculation any further, as we are only 
interested in the qualitative behavior of Pf . We just note that these factors can be organized 
in a scheme resembling Pascal's triangle. 

Next we give the proof of the approximation in eq. ( 130|) . We investigate the integral 
corresponding to the level lj for which the energy difference to the initial level is maximal. 
For convenience we define e := |Aej.j| = |a'(a;)| (a(x) is defined in eq. (EH}) and /i := lA/^. 

In the present case we assume that the spectrum is peaked at a sufficiently high frequency 
Qq, with the property Qq » e. In addition the width of the peak is assumed much smaller 
than Q . In this case Q e [0, e] hardly makes any contribution and can be ignored for finding 
an upper bound to \I n (E,t)\. 

In the following we set the lower limit of integration to the moment to when the field 
was switched on. Because the derivative of the field up to any order is zero at t < to it is 
clear that the lower boundary does not make a contribution and sin n (nt + 0) is replaced 
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by sin(fit + <p)/Q n and similar for the modified cosine. From eq. ( 1B8j) we obtain the upper 
bound 



\In(E,t)\ < - 
7T 



dtt- 



(C3) 



which is equal to the maximally attained value of the oscillating function I n (E,t). Noting 
that the peak of S(Q) is much larger than its halfwidth we obtain the order estimate 



\i n (E,t)\ = -o(n n ) , 

7T 



(C4) 



which we use in the proof below. 

Now we show that e << f2 indeed implies the approximation made in eq. f )30|) for n = 1 
and extend to n > 1 later. The approximation clearly is valid whenever 

»ti 



dxa'(x)B(x) 



« \a(x)B(x) 



(C5) 



holds; the abbreviations are defined in eq. ( |29|) . In order to motivate the strategy of the proof 
we note that multiple integrals of b(x) are oscillating functions, that is, strictly the lower 
bound of their magnitude is zero. Therefore we cannot proceed to construct an upper/lower 
bound for the left/right side of eq. (105|) in order to prove it. We then use the estimate 
eq. (lC4p . noting that the inequality eq. ( 105j) becomes invalid in small intervals around 
isolated zeros of B(x). However, these intervals occupy a fraction of O(e/Qo) of the whole 
integration time and noting that within these intervals we have no incremented by 1, which 
multiplies an extra factor e/Qo to the wavefunction, the contributions from these intervals to 
the population are by a factor (9[(e/f2 ) 2 ] << 1 smaller than from the remaining intervals. 
Using eq. (104[) we immediately get 



dxa'(x)B(x) 



A 

7T 



- oi m 



(C6) 



Next we have to obtain a suitable lower bound to the right hand side of eq. (IC5j) . We note 
that introducing z(x) := — exp(— %\i f Q dx'E(x')) we obtain dz(x)/dx = inE(x)z(x) = b(x). 
B(x) can then be integrated analytically to give 



B(x) = - M 



exp 



dt'E(t') 



(C7) 



The exponent evaluates to j*dt'E(t') = 0{A/nVL Q ) w O(e/Q ) 

« 1 when we consider field 

strengths of the order of the one from the dipole resonance condition (see IIIIBTT) . Due to 
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the fact that the length of 1 — exp(— i<f>) increases monotonically with increasing £ [0, it] 
we are safe to postulate 



1 — exp 



px "1 rX 

ifi I dt'E(t') « fi dt'E(t' 
Jo . Jo 



(^) . 

7T 



(C8) 



Inserting this result in eq. (1CT[) we get 



.4 



\b(x)\ = -o fa 1 ) , 

7T 

which tells us that |/b(1,^)| ~ l-M-E^)! is a good approximation in our case. 
Putting it all together we immediately arrive at the desired result, 

e << n n ■ 



(C9) 



(CIO) 



In case of n > 1 we note the relation \I B (k,t)\ < \Ik(E,t)\ for any k, with decreasing 
difference \Ik(E,t) \ — \I B (k,t)\ for increasing Q . We again use the estimate from eq. (jC4p 
for both sides of the inequality to arrive at 

,n -l 



c'K) 



Si 



<< 



li 



(Gil) 



which clearly is equivalent to eq. flClOj) and therefore the approximation in eq. f l30]) holds 
for general Uq. 
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Figure captions 

Figure 1. (Color online) Maximum population P f max (<i) in the final level of a 2LS under 
constant field, plotted against E/A (upper panel) and against the dipolar detuning (lower 
panel). The levels are denoted i (initial) and f (final). The resonance field strength is set 
to Aq = Atfi/A/ifi = 1, and two values of the ratio /ifi/A/ig are used as indicated. The 
maximum detuning is given by A/xg/2/ifi, leading to the (visible) right-hand cutoffs in the 
lower graph. 
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FIG. 1: 
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